Directed flow in ultrarelativistic heavy-ion collisions 
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We study the generation of the directed flow in the hydrodynamic expansion of hot matter formed 
in ultrarelativistic heavy-ion collisions at y/s — 200GeV. The experimentally observed negative 
directed flow in a wide range of central pseudorapidities is reproduced, assuming that the fireball is 
tilted away from the collision axis. The tilt of the source is consistent with a preferential emission 
in the forward/backward hemisphere from forward/backward participating nucleons. The model 
reproduces the experimentally observed scaling of the directed flow when going from Au-Au to 
Cu-Cu systems. 
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I. INTRODUCTION 

The appearance of the transverse, azimuthally asym- 
metric flow is one of the key observations in the physics 
of relativistic heavy-ions [l[ . It proves that a collectively 
expanding, strongly interacting medium is formed in the 
course of the reaction. A number of observables have 
been studied, both in experiments, and in model calcu- 
lations, in order to unravel the properties of the dense 
hot matter created in the collisions. The production of 
particles with soft momenta can be interpreted as a ther- 
mal emission of particles from fluid elements moving with 
some collective velocity field Q. Relativistic hydrody- 
namics describes quantitatively the development of the 
collective velocity from pressure gradients in the fireball 

For non-central collisions the interaction region is az- 
imuthally asymmetric and, as a result of the collective 
expansion of matter, azimuthally asymmetric emission 
of particles takes place. The effect can be quantified in 
terms of Fourier coefficients in the expansion of the mea- 
sured particle spectra 
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The elliptic flow coefficient V2 is known to be a very sen- 
sitive probe of the pressure in the system at early stages 

Directed flow, quantified by the coefficient v\ is also 
measured at energies available at the BNL Relativistic 
Heavy Ion Collider (RHIC). The coefficient V\ of the di- 
rected flow is zero at zero rapidity for collisions of sym- 
metric nuclei, but it increases when moving to forward or 
backward pseudorapidities r\. Its size and sign has been 
the subject of many studies at lower energies, where it 
is dominated by nucleon flow Q. At RHIC energies the 
spectator nucleons have positive directed flow (y\ > 
for 77 > 0), resulting from the deflection of the spectators 
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during the collision. On the other hand, matter in the 
fireball shows a significant negative (anti-flow) compo- 
nent for pseudorapidities — 4 < rj < 4 and for centralities 
c = - 80%, both for Cu-Cu and Au-Au collisions [1- 
[Toj . A striking characteristic of the measured directed 
flow is the large negative value of V\ even at the high- 
est energy y/s = 200GcV. Another important observa- 
tion is the scaling of the measured directed flow with the 
size of the system. The coefficient v\ is the same for 
both systems (Au-Au or Cu-Cu), for the same centrality. 
Whereas, a scaling of i>2 with the density of the fireball 
has been observed [TlT | . 

Transport models of nuclear reactions, describe the 
directed flow at lower energies, but generally underpre- 
dict the amount of anti-flow at RHIC energies at central 
rapidities 12]. On the other hand, these calculations 
predict a large negative flow at large rapidities. Some 
calculations 13] predict significant anti-flow around cen- 
tral rapidities, larger than observed, but yield a positive 
flow for I77I > 3, unlike observed in experiments. It has 
been noticed that the appearance of negative directed 
flow around central rapidities could be an effect of the 
softening of the equation of state [13] • This effect is 
called the third flow component. Hydrodynamic calcula- 
tions incorporating such effects yield a negative elliptic 
flow at central rapidities and a positive directed flow at 
larger pseudorapidities, unlike the experimental data. A 
hydrodynamic calculation with initial conditions from a 
microscopic model gives the correct sign, but a too small 
magnitude of the directed flow for central rapidities (15| . 

There are two effects leading to a negative directed 
flow in the models. The first one is the shadowing of the 
fireball matter by the spectators, which can give a sub- 
stantial negative v\ for \rj\ > 4. The second one is the 
build up of the flow away from the collision axis due to 
a tilt of the source. The description of the magnitude 
of the directed flow at different centralities for central 
rapidities requires, a tilt of the source of the right magni- 
tude as function of the impact parameter, and a sufficient 
amount of collectivity to generate the flow. We consider 
two initial conditions for the hot source in the hydro- 
dynamic evolution. The first one is the most commonly 
used initial conditions in 3+1-dimensional (3+ ID) calcu- 
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lations incorporating a Bjorken flow in the longitudinal 
direction and a shift in space-time rapidity due to the 
local imbalance of the momentum The second one 
assumes that the initial density results from a superpo- 
sition of the energy density radiated by the color sources 
in the target and the projectile. The preferred emission 
in the pseudorapidity hemisphere of the emitting charge 
results in a tilt of the source for non-central collisions. 
We show that the second choice leads to a satisfactory 
description of the directed flow generated in heavy ion 
collisions at yfs — 200GeV for central rapidities. 



II. INITIAL CONDITIONS AND EARLY FLOW 

Hydrodynamic evolution in 3 + ID at RHIC ener- 
gies is performed in the proper time r = \/t 2 — z 2 . 
The densities are defined as functions of the trans- 
verse plane (x-y) coordinates and the space-time rapidity 
m = h log ((t + z)/(t — z)). The hydrodynamic model 
requires some initial density and initial flow profile to be 
chosen at the initial time tq. Although some guidance 
from microscopic models of elementary collisions is pos- 
sible in the choice of the initial conditions, there is still a 
vast choice of initial energy density profiles that are used 
in the simulations. Also, most of the calculations assume 
a Bjorken initial flow in the longitudinal direction and no 
initial transverse flow 



u^(To,x,y,r)\\) = (cosh 77||, 0,0, stimuli ) 
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In the hydrodynamic evolution the generation of 
the left-right asymmetry of the flow in the two (for- 
ward/backward) halves of the reaction plane requires the 
presence of an asymmetry in the initial distributions. For 
non-central collisions the azimuthal asymmetry of the in- 
teraction region results in a non-zero initial eccentricity 
in the transverse plane that gives rise to the collective 
elliptic flow. For collisions of symmetric nuclei and ne- 
glecting the fluctuations, the odd components of the de- 
composition in Fourier coefficients vanish at space-time 
rapidity zero. At forward and backward rapidities an im- 
balance between the contributions from the target and 
the projectile to the initial source can result in a left- 
right deformation of the source in the transverse plane 
and/or in an asymmetric initial flow, that could generate 
collective directed flow. 

The observation of non-zero directed flow implies, that 
the symmetry in the reaction plane is indeed broken, ei- 
ther in the initial flow or in the initial density, or both. 
We tried to reproduce the observed directed flow assum- 
ing an asymmetric initial flow, different from the Bjorken 
one @, but without success. Therefore in the following 
we assume a Bjorken initial flow and study the effect of 
asymmetric initial densities on the directed flow. 

Starting with a factorized initial energy density in the 
transverse plane and in the space-time rapidity 



the subsequent evolution remains symmetric with respect 
to the m axis in the reaction plane (n\-x), with the conse- 
quence that the directed flow is exactly zero. Factorized 
initial conditions in the Glauber Model imply that all 
the participant nucleons or binary collisions contribute 
in a similar way to the total density, with a longitudinal 
profile proportional to /foil). 

Modifications of the symmetric distribution ([3]) could 
happen during the formation of the initial thermalized 
state. Momentum imbalance between left and right go- 
ing participants at a given point in the transverse plane 
results in a nonzero total momentum of the matter. The 
longitudinal distribution is shifted in space-time rapidity 
[H EH by the value of the center of mass rapidity of the 
fluid 

1 { N + +N-+v N (N + -N_) \ 
Vsh = 2 l0g { N+ +N--v N (N + -N-) ) ' (4) 

where N + and 7V_ are the densities of participants from 
the two nuclei and Ujy is the velocity of the incident nu- 
clei. 

N + (x, y) = T(x - b/2, y) (f - exp(- ^±^ ) 

7V_ (x, y) = T(x + b/2, y) (l - exp(- ^ '^ v) ) }f, , 
where a is the cross section, 



T(x,y) = / dzp(x,y,z) 



(6) 

is the thickness function calculated from the Woods- 
Saxon density of colliding nuclei 

Po 



p(x,y,z) 



I + exp \ {\J x 2 + y 2 + z 2 - R A )/a 



(7) 



t(v\\,x,y) = p(x,y)f{v\\) 



(3) 



We use the same parameters as in Ref. [Bj], where a sat- 
isfactory description of spectra and femtoscopy data for 
Au-Au collisions at yfs = 200GeV has been obtained. 
The initial energy density takes the form 0, EH 

e(r ) - eo/fau - Vsh) [{N+(x, y) + N-{x, y)) (I - a) 
+ 2aN bin (x,y)]/N . (8) 

The relative contribution from binary collisions is a — 
0.145, with N bin (x, y) = aT(x - b/2, y)T(x + b/2, y). In 
the following we call these initial conditions shifted initial 
conditions (Fig. [T]). 

The form of the initial energy density profile /foil) 
is adjusted to reproduce the measured charged parti- 
cle distribution in pseudorapidity. The resulting width 
of the initial distribution depends on the chosen equa- 
tion of state, initial time tq and shear viscosity [5l Il7|. 
In the following we use ideal fluid hydrodynamics with 
To = 0.25fm/c and a realistic, hard equation of state [l8j . 
which requires 

/foil) = ^p (- ^^J° )2 e(\ m \- r, )) (9) 
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FIG. 1. (Color online) Contour plot of the initial pres- 
sure p{r), x,y — 0) in the fireball for the shifted densities 
(Eq. [3J. Solid lines correspond to the pressure of 9, 3 and 
fGeV/fm 3 for Au-Au collisions (impact parameter b = llfm) 
and dashed lines to the pressure of 3 and f GeV/fm 3 for Cu- 
Cu collisions (6 = 7.6fm). The arrows represent the gradient 
(—d v p/To, —d x p) for Au-Au collisions in arbitrary units. 
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FIG. 2. (Color online) Initial profile in the longitudinal 
(space-time rapidity) direction. The symmetric function 
/(77||) is composed from two contributions /+ and /_ rep- 
resenting the emission from forward and backward going par- 
ticipant nucleons. 



with a plateau of width 2rjo — 2.0 and a n = 1.3. 

A different type of initial conditions studied in this 
work assumes a preferred emission from participating nu- 
cleons in the same hemisphere. Instead of a symmetric 
distribution of matter in space-time rapidity given by the 
function f(r]n) in Eq. ([3]), we assume that the deposited 
energy depends on the rapidity of the emitting participat- 
ing nucleon. Such a distribution depending on the rapid- 
ity difference between the emitting charge and the emit- 



FIG. 3. (Color online) Same as in Fig. [IJbut for tilted initial 
conditions (Eq. [T3)) . 



ted gluon is assumed in some phenomenological models 
19]. However, there is no direct measurement of the con- 
tribution to soft particle production from a single forward 
or backward moving charge. A phenomenological anal- 
ysis is possible, by comparing multiplicity distributions 
in pseudorapidity for different asymmetric systems or by 
studying multiplicity correlations in different pseudora- 
pidity intervals. These studies indicate that a preferred 
emission for rapidities close to the rapidity of the partic- 
ipating charge occurs [20l - [23| . In the wounded nucleon 
model of nuclear collisions, such correlations can be un- 
derstood as due to a specific distribution of soft particles 
produced by each participant nucleon. Nucleons from the 
projectile (with positive rapidity ys = ln(\/s/mjv) > 0) 
emit more particles in the forward (77 > 0) than in the 
backward hemisphere. The form of the extracted charged 
particle distribution can be approximated by the function 



Mr?) 



V + Vn 



(10) 



in the interval [— r) m , r) m ], where r\ m = yb — r/ a defines 
the range of rapidity correlations, at y/s = 200GeV it is 
rj m ~ yb — 2 ~ 3.36. The origin of the shift in rapidity 
ij s ~ 2 is not understood [2(J Hlj . For practical purposes 
we can treat it as a phenomenological parameter. Parti- 
cle production in the remaining pseudorapidity intervals 
close to the fragmentation regions [rj m , Ub] and [—yb, —Vm] 
cannot be reliably described in a hydrodynamic model 
anyway. Within the framework of relativistic hydrody- 
namics, we are interested in describing the main charac- 
teristics of the soft part of particle spectra in the central 
region —3.5 < r\ < 3.5 and, in particular, the directed 
flow. There is another reason why the phenomenological 
estimates of the emission of particles from participant 
nucleons [2(| Hl[ cannot be directly translated into the 
initial conditions for hydrodynamics, that we are inter- 
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ested in. Refs. [20, |2l[ study particle distributions and 
correlations in the final state, whereas, we know that in 
realistic hydrodynamic simulations the matter distribu- 
tion in space-time rapidity evolves during the expansion 
of the fireball p| [I?} , also statistical emission broadens 
the distribution in pseudorapidity. ft means that the ini- 
tial profile f(f]u) is significantly narrower than the final 
charged particle distribution The correlation func- 
tions in pseudorapidity 2l|] can be modified due to the 
longitudinal transport and the generation of transverse 
flow as well. We propose as a phenomenological ansatz, 
inspired by the observations in Refs. [20l - l23l |. that the 
initial energy density of matter produced by a single par- 
ticipant nucleon of rapidity yi> is proportional to 



f+(v\\) = f(V\\)fF(V\\ 



(11) 



where f(r)\\) is the initial longitudinal profile © fitted to 



reproduce and 



fp(ri\\ 





1 



-Vm < V\\ < Vn 
Vm < V\\ 



(12) 



The initial energy density of the fireball is constructed 
as a sum of three terms originating from the forward or 
backward moving participant nucleons and from the bi- 
nary collisions that are assumed to contribute in a sym- 
metric way 

e(r ) = eo [2 (JV+(x, y)f+(r)\\) + N-(x, y)/-(r?||)) (1 - a) 
+ 2aN bm (x,y)f( m )] /N . (13) 

The net result of the difference between forward and 
backward emission is a tilt of the source in the x-r]\\ plane 
(Fig. This breaks the symmetry in the longitudinal 
direction and generates nonzero directed flow in the ex- 
pansion. 

Hydrodynamic equations in 3 + ID 







(14) 



constitute four independent equations, which, together 
with the equation of state, determine the evolution of 
the energy density e, of the pressure p, and of three inde- 
pendent components on the fluid velocity. The fluid four 
velocity can be written in the form 



(7 cosh Y, u x , u y , 7 sinh Y ) 



(15) 



u x and u v are the components of the transverse veloc- 
ity, 7 = J\ 



Uy and Y 



l-v 2 



is the fluid 



rapidity. The densities are functions of the proper time 
r, the space-time rapidity r)\\ and the transverse coordi- 
nates x, y. The equations in the expanded form can be 
found in Q. At early times, the velocities on the right 
hand side of the equations can be approximated by the 



initial velocities u x = 0, u y =0, Y = m. The two ac- 
celeration equations for the velocity components in the 
reaction plane take the form 



d T u x 



1 



e + p 



-d x p 



d T Y 



r(e + p) 



d m P 



(16) 



In 3 + ID hydrodynamic evolution, the lack of Bjorkcn 
invariance results in a nonzero longitudinal acceleration. 
The fluid rapidity Y becomes larger than the space-time 
rapidity m. In Figs. [1] and [3] are shown the vector fields 
of the initial pressure gradients (-^-dr^p, d x p). For the 
shifted initial conditions the gradient in the central region 
of the fireball is in the transverse direction, and mainly 
transverse flow is generated in the early stage. This is 
due to the existence of an approximate Bjorken plateau 
for central rapidities in the initial stage. The situation is 
different for tilted initial conditions (Fig. the acceler- 
ation in the tilted source is anti-correlated in the trans- 
verse x and longitudinal 7711 directions. The matter that 
is accelerated to positive rapidities is preferably acceler- 
ated in the negative x direction. In the same Figures are 
shown the lines of constant pressure for the initial fire- 
ball created in Cu-Cu collisions for the same centrality. 
The deformation for the shifted fireball (Fig. [T]), or the 
tilt for the tilted fireball (Fig. [3]), in the smaller system 
is very similar as in the larger system. This generates a 
similar directed flow in the two systems irrespective of 
their sizes. 



III. RESULTS 

The hydrodynamic equations are solved numerically 
for the two sets of initial conditions, the shifted ini- 
tial distributions (jHJ and the tilted initial conditions 
(fT3"|). The parameter eo is chosen to reproduce particle 
spectra and multiplicities in central collisions. We use 
e = 107GeV/fur 3 and e = 65GeV/fm" 3 for Au-Au 
and Cu-Cu collision respectively, and a freeze-out tem- 
perature of Tp — 150MeV. This gives a satisfactory de- 
scription of the spectra in collisions up to centralities of 
50% @. 

The distribution of charged particles in pseudorapid- 
ity is shown in Fig. |U The results obtained from the 
two initial conditions are almost indistinguishable on the 
plot. Both initial conditions lead to similar results for 
transverse momentum spectra of particles, interferome- 
try radii, and elliptic flow, as well. The comparison with 
experimental data can be found in Ref. [5[, giving sat- 
isfactory results. The calculated elliptic flow overshoots 
the experimental data since we do not take viscosity ef- 
fects into account 26]. The parameters of the initial 
profile in space-time rapidity /(?7||) are adjusted to repro- 
duce, as closely as possible, the experimental results on 
pseudorapidity distributions, and transverse momentum 
spectra at non-zero rapidities Changing the width 
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FIG. 4. (Color online) Pseudorapidity distribution of charged 
particles for centrality classes — 6%, 15 — 25%, calculated 
for the shifted and tilted initial conditions (solid and dashed 
lines respectively) compared to PHOBOS Collab. data (dots) 
[23 |. The squares represent the BRAHMS Collab. data for 
centrality - 5% [Hf. 
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FIG. 5. (Color online) Average flow in the x direction as 
function of space-time rapidity for different evolution times 
for the Hirano-Tsuda shifted initial densities (Eq. [8}. 



r/Q of the plateau in the initial profile /(?7||), within a 
range compatible with the observed pseudorapidity dis- 
tributions, does not change the results for the directed 
flow. Only taking the unrealistic value 770 = for shifted 
initial conditions, causes the wiggle of positive v\ disap- 
pear at central pseudorapidities. 

By breaking the symmetry in the longitudinal direc- 
tion, some directed flow can be generated. In Figs. [5] 
and |B] is shown the development of the asymmetric flow 
in the x direction at different times. The average veloc- 
ity in the x direction is calculated for a given space-time 
rapidity and time 
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FIG. 6. (Color online) Same as in Fig. [5] but for tilted initial 
conditions (Eq. 113}) . 



As mentioned above, during the evolution the Bjorken 
flow Y = is modified, but the final flow conserves a 
strong correlation between Y and r\\\ |f|. Therefore the 
directed flow as functions of space-time rapidity r}\\ as 
shown in Figs. [5] and [6] reflects qualitatively the final 
flow of particles. 

The first observation is that the directed flow is built 
up for an extended time during the expansion, in the 
first 3fm/c. There is a noticeable difference between the 
evolution from the shifted and the tilted initial condi- 
tions. For the shifted initial conditions (Fig. [5]) in the 
central space-time rapidity region a positive directed flow 
develops gradually in the first lfm/c. For large rapidi- 
ties a negative flow is generated in the first 3fm/c. The 
directed flow changes sign around i] \\ — 2, a behavior 
very different from the one observed in the experiment. 
The appearance of the wiggle in the dependence of v\ on 
pseudorapidity is a consequence of the form of the initial 
profile of the pressure (Fig. [1]). The gradient of the pres- 
sure has a significant deflection from the transverse direc- 
tion only in the forward/backward rapidity regions. On 
the other hand, the tilted initial condition has a smooth 
tilted gradient of the pressure that gives a negative di- 
rected flow in a broad range of space-time rapidities (Fig. 
[5]). The anti-flow increase for the first 3fm/c. The neg- 
ative directed flow of the fluid, increasing with rapidity, 
leads to a similar pattern in the directed flow of the final 
particles. 

The hydrodynamic evolution is continued until freeze- 
out, where particle emission from the freeze-out hyper- 
surface takes place. Statistical emission and resonance 
decays are performed using the event generator THER- 
MINATOR H3j. In Figs. [7] and M are shown the results 
for three representative centralities, for which experimen- 
tal data for Au-Au collisions have been published. Other 
experimental data show that in Cu-Cu interactions al- 
most the same directed flow is generated as in the larger 
system if the centrality is chosen to be the same (luj . 
The directed flow for the shifted initial conditions shown 
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FIG. 7. (Color online) Directed flow at different centralities 
from hydrodynamic calculations with shifted initial conditions 
(Eq. [8]) for Au-Au and Cu-Cu collisions (solid and dashed 
lines). Experimental data are from the PHOBOS Collabora- 
tion for c = - 40%, and from the STAR Collaboration 
fiol ] for the three other centrality classes. 



FIG. 8. (Color online) Directed flow in Au-Au (thick solid 
lines) and Cu-Cu (dashed lines) collisions at different cen- 
tralities from tilted initial conditions (Eq. I13|) . compared to 
experiment [13] . The shaded band between the thin and 
thick lines represents the increase of the magnitude of the 
flow if Eq. p8p is used for the initial density, including also 
asymmetric contributions from binary collisions. 



in Fig. [7] has an incorrect dependence on pseudorapid- 
ity. For all the three centrality classes the flow is positive 
in the central rapidity region, and switches to anti-flow 
in the very forward and backward regions. On the same 
plots are shown the results for Cu-Cu collisions in similar 
centrality classes. The flow is similar as for the larger sys- 
tem, and hence different than observed experimentally. 
The similarity between the flow in Au-Au and Cu-Cu 
systems reflects the similarity in the initial density pro- 
files (Fig. Q). 

In Fig. [8] are shown results for the directed flow of 
charged particles emitted after a hydrodynamic evolu- 
tion from the tilted initial fireball (Eq. [T3|) . The thick 
solid lines represent the results for Au-Au collisions. The 
experimental data are reproduced in the central rapidity 
region in a satisfactory way. For the semi-peripheral and 
peripheral collisions the large anti-flow at large pseudora- 
pidities is not reproduced by the model. This kinematic 
region is at the limit of applicability of the hydrodynamic 
model, assuming the collective expansion of a thermal- 



ized fluid. The directed flow in the fragmentation region 
can have a different origin, like in the transport models 
12]. On the panels in Fig. [5] are shown the results for 
the Cu-Cu system in the same centrality classes. Again, 
as expected from the similarity of the tilt of the source 
in the two systems (Fig. [3]), the final directed flows for 
the two systems almost overlap. To test some of the 
uncertainty in the choice of the initial conditions for the 
hydrodynamic evolution, we take as the initial density an 
ansatz where also the contribution from binary collisions 
is asymmetric 



e(r ) = 2e 



N + U(m)+N_f_(m) 



N- 

[(l-a)(N++N-) + 2aN bin ]/N . (18) 

The tilt is stronger in this case and the magnitude of 
the directed flow bigger (thin solid lines in Fig. [5]). The 
shaded bands in the figure represent the uncertainty of 
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FIG. 9. (Color online) Directed flow coefficient as function of 
pseudorapidity for tilted initial conditions (Eq. I13[) in Au-Au 
collisions. The solid and dashed lines represent the calcula- 
tions with initial times to = 0.25 and lfm/c respectively and 
with an initial profile fp with rj m = yj, — 2, the dotted line is 
for to — 0.25fm/c but rim = yb, and the dashed-dotted line 
represents the result of a calculation where the initial entropy 
density is proportional to the density of part icipants from the 
Glauber Model. Data are from Refs. [si. TTol| . 



the model related to this assumption. This uncertainty 
is the largest among those that we tested. 

In the following we study the effect of other details 
of the model on the final directed flow. First, a differ- 
ent initial time for the hydrodynamic evolution is cho- 
sen To = lfm/c. The initial time is given by the rate 
of the initial thermalization processes, that lead to the 
formation of a dense, almost perfect fluid. The precise 
mechanisms of thermalization are not known and it is 
instructive to test the influence of the value of the ther- 
malization time on the build up of the directed flow. In 
Fig. [9]is shown the result for the centrality class 5 — 40%. 
The directed flow is reduced when compared to the pre- 
vious calculation using tq = 0.25fm/c. The fluid that 
starts to expand at tq = lfm/c has a smaller pressure. 
In fact, the experimental points for central rapidities lie 
between the two lines representing the two calculations 
starting at different thermalization times. 

Another effect worth to be studied is the dependence 
on the details of the form-factor f±(rjn), describing the 
emission from a single participant nucleon. The range 
of forward-backward correlations is set by the parameter 
r] m . We make another calculation using rj m = yb. With 
this choice, the emission from a forward going participant 
nucleon decreases when going from rapidity yj, to zero at 
the rapidity —yb- The initial tilt of the source is smaller 
than for the choice rj m = 2/6 — 2, used before. The effect of 
this change on the directed flow is small and the results 
are close to the experimental data at central rapidities 
(dotted line in Fig. 

A third effect, that can be tested, is the assumption 
that the initial energy density is proportional to the den- 



sity of wounded nucleons (binary collisions) as in Eq. 
([TBI . Assuming that this proportionality applies to the 
entropy density s instead, we have 

s(ro) = so [{N + (x,y)U( m ) + N_(x,y)f-( m )) (1 - a) 

+ 2aN Un {x,y)f( m )] /N . (19) 

The directed flow is very similar to the one obtained us- 
ing the energy density initial profile (fT3)) (dashed-dotted 
line in Fig. |9]). We notice that some of the details of the 
initial conditions can influence the final directed flow. 
This situation resembles the conclusion from the ellip- 
tic flow studies, where a crucial ingredient is the initial 
eccentricity of the fireball. 

IV. SUMMARY 

We study the formation of the directed flow in the 
3 + ID hydrodynamic expansion of the fireball created in 
heavy-ion collisions at the highest energies at RHIC. The 
directed flow of charged particles has been measured as 
function of pseudorapidity, finding a substantial negative 
flow We use two different initial conditions for 

the evolution. The first one (Eq. [8| is quite commonly 
used in hydrodynamic model calculations [3]. It incor- 
porates a shift of the densities in the initial fireball, due 
to the local imbalance of the longitudinal momentum. 
Initial densities of the second type are constructed as a 
sum of contributions from forward and backward going 
participants (Eq. I13[) . The asymmetry in the emission 
from individual participants leads to a tilt of the source. 
Model calculations, incorporating a hydrodynamic ex- 
pansion stage, particle emission at freeze-out, and res- 
onance decays, indicate that the second type of initial 
conditions can reproduce the sign and the magnitude of 
the observed directed flow at central rapidities. 

The large negative flow close to the fragmentation re- 
gions is of different origin and cannot be described in 
our calculation. Deviations from the Bjorken flow in the 
initial conditions [28[ could influence the results. In par- 
ticular, they could lead to a different flow pattern for 
baryons and for the bulk of the matter, if baryons do not 
follow the Bjorken flow. The directed flow is expected to 
decrease with the collision energy for two reasons. With 
increasing kinematic range of rapidities the tilt of the 
source goes down, also the contribution of binary colli- 
sions to the Glauber Model density is believed to increase 
with the energy, with a similar consequence. Finally let 
us note, that our calculation reproduces the experimen- 
tally observed similarity of the flow in Au-Au and Cu-Cu 
collisions at the same centrality. 
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